Synthesis and characterization of gold(I) thiolate derivatives and bimetallic complexes for HIV inhibition

Introduction: The human immunodeficiency virus (HIV) remains a significant global health concern, with a reported high infection rate of 38.4 million cases globally; an estimated 2 million new infections and approximately 700,000 HIV/AIDS-related deaths were reported in 2021. Despite the advent of anti-retroviral therapy (ART), HIV/AIDS persists as a chronic disease. To combat this, several studies focus on developing inhibitors targeting various stages of the HIV infection cycle, including HIV-1 protease. This study aims to synthesize and characterize novel glyco diphenylphosphino metal complexes with potential HIV inhibitory properties. Method: A series of new gold(I) thiolate derivatives and three bimetallic complexes, incorporating amino phosphines and thiocarbohydrate as auxiliary ligands, were synthesized using procedures described by Jiang, et al. (2009) and Coetzee et al. (2007). Structural elucidation and purity assessment of the synthesized compounds (1–11) were conducted using micro-analysis, NMR, and infrared spectrometry. Results and Discussion: Using molecular modeling techniques, three of the metal complexes were identified as potential HIV protease inhibitors, exhibiting strong binding affinity interactions with binding pocket residues. These inhibitors demonstrated an ability to inhibit the flexibility of the flap regions of the HIV protease, similar to the known HIV protease inhibitor, darunavir. This study sheds light on the promising avenues for the development of novel therapeutic agents against HIV/AIDS.


Introduction
Communicable diseases represent significant global challenges, impacting human security, exacerbating poverty, and fostering economic instability (Beagleholeet al., 2011).Among these diseases, the human immunodeficiency virus (HIV) is reported to have a high infection rate, with a global prevalence of 38.4 million in 2021, and there were an estimated 2 million new infections and approximately 700 000 HIV/AIDS-related deaths in the same year (UNAIDS and Sheet, 2021).Despite advances in treatment, HIV/ AIDS remains a chronic condition, necessitating ongoing efforts to develop novel therapeutic approaches (Levitt et al., 2011).
HIV protease is among the biomolecules targeted in the design of HIV therapeutics (Arrigoni et al., 2023).It is involved in the cleavage of viral polyproteins, leading to viral protein maturation (Chuntakaruk et al., 2024).Therapeutically blocking HIV protease therefore halts viral maturation, thus presenting HIV protease as a prominent target for HIV drug discovery efforts.Several HIV protease inhibitors have been developed over the years (King et al., 2004;Ghosh et al., 2006;Chuntakaruk et al., 2024).Among these is darunavir, whose binding is characterized by strong interactions between the 3(R),3a(S),6a(R)-bistetrahydrofuranyl (bis-THF) moiety of darunavir and the main-chain atoms of aspartates 29 and 30 (King et al., 2004;Ghosh et al., 2006;Chuntakaruk et al., 2024).
Advancements in metal-based antiviral therapeutics, particularly involving gold (Au) and silver (Ag), have shown promise in enhancing safety and minimizing side effects (King et al., 2004;Ghosh et al., 2006).Metal-drug synergism, where drug stability is enhanced by coordination with metal ions, offers advantages such as improved pharmacokinetics and drug metabolism (Pettersen et al., 2004;Trott and Olson, 2010;ChemAxon, 2013;Kožíšek et al., 2014).Gold-based therapeutics, in particular, have demonstrated efficacy in inhibiting viral entry by interacting with Gp120 and CD4 (Hanwell et al., 2012).
Despite uncertainties surrounding the precise mechanisms of action, studies have shown that gold compounds effectively block viral entry and inhibit HIV proteases (Hanwell et al., 2012;Le Grand et al., 2013).Notably, gold(I) phosphine compounds have exhibited inhibitory activity against HIV-1 reverse transcriptase (RT) and protease (PR) without causing toxicity to immune cells (Le Grand et al., 2013).Metal-carbohydrate compounds have also shown potential in reducing toxicity and enhancing biocompatibility, offering promising avenues for therapeutic development (Ramharack et al., 2023).However, concerns regarding side effects, such as skin discoloration upon light exposure, have been reported with gold-based therapeutics (Agoni et al., 2018a).
Although antiretroviral therapy (ART) has extended the lifespan of HIV patients, there remains a need for therapeutics capable of effectively inhibiting HIV-infected cells (Olotu et al., 2019).To address this need, pharmaceutical companies are increasingly leveraging computer-aided techniques to predict drug candidates, reducing research costs and failures attributed to poor pharmacokinetic properties (Agoni et al., 2022).In silico molecular docking, molecular simulation, and pharmacokinetic analyses offer valuable insights into ligand-target protein interactions, aiding in the identification of potential therapeutic agents (Salomon-Ferrer et al., 2013;Maier et al., 2015).
In this study, we seek to synthesize and characterize novel glyco diphenylphosphino metal complexes with the potential to combat HIV.Leveraging advanced molecular modeling techniques, we aim to elucidate the mechanisms underlying the anti-HIV activities of these compounds at a molecular level.By unraveling these intricate interactions, we seek to contribute to the development of more effective therapeutic strategies against HIV.

Instrumentation
All NMR spectra of CD 3 OD and CDCl 3 were recorded using a Bruker UltraShield Spectrometer (400 MHz) at room temperature.The 1 H and 13 C{ 1 H} chemical shifts were referenced to the residual signals of the protons or carbons of the NMR solvents and are quoted in ppm, i.e., CD 3 OD at 4.79 and 49.00 ppm for 1 H NMR and 13 C{ 1 H} spectra, respectively, and CDCl 3 at 7.24 and 77.00 ppm for 1 H and 13 C{ 1 H} spectra, respectively.The infrared spectrum was recorded using a Bruker Tensor 27 Spectrometer fitted with an ATP-IR probe and PerkinElmer FTIR Spectrum BX.Elemental analysis was performed using a Vario Elementar III Microcube CHNS Analyzer at the University of Johannesburg, South Africa.

Synthesis of [2-(diphenylphosphino)ethyl] amine gold(I) chloride (1)
Title compound 1 was prepared according to a previous procedure (Jiang, et al., 2009) from [Au(tht)Cl], which was obtained by reacting H [AuCl 4 ].3H 2 O with tetrahydrothiophene (tht) in ethanol.A dichloromethane solution of [Au(tht)Cl] (0.26 g, 0.87 mmol) was added to a solution of 2-(diphenylphosphino)ethanamine (L1) (0.20 g, 0.87 mmol) in dichloromethane (10 mL) and stirred for 3 h.The solvent was reduced to half, and hexane was added until a white precipitate was formed.The white solid obtained was filtered and dried in vacuo.

Synthesis of [acetylated-2-gluconamido ethyl (diphenylphosphino)] gold(I) chloride (5)
To a solution of complex 3 (0.12 g, 0.18 mmol) and acetic anhydride (0.10 mL, 1.10 mmol) (1:1 ratio to the hydroxyl group) in pyridine (2 mL), a catalytic amount of DMAP was added.The solution was stirred at room temperature for 18 h.The mixture was then diluted with dichloromethane (20 mL) and washed with 1M HCl (aq) solution (5 × 30 mL).The organic layer was dried using anhydrous MgSO 4 , after which the solvent was pumped off using a rotary evaporator to obtain a white solid.Yield: 0.16 g (84%).

Synthesis of [acetylated-3-gluconamido propyl (diphenylphosphino)] gold(I) chloride (6)
To a solution of complex 4 (0.10 g, 0.15 mmol) and acetic anhydride (0.08 mL, 0.75 mmol) (1:1 ratio to the hydroxyl group) in pyridine (2 mL), a catalytic amount of DMAP was added.The solution was stirred at room temperature for 18 h.The mixture was then diluted with dichloromethane (20 mL) and washed with 1M HCl (aq) solution (5 × 30 mL).The organic layer was dried using anhydrous MgSO 4 , after which the solvent was pumped off using a rotary evaporator to obtain a white solid.Yield: 0.11 g (85%).

Molecular docking
AutoDock Vina (Trott and Olson, 2010) was used for the molecular docking of the synthesized mono-and bimetallic glyco diphenylphosphino gold(I), palladium(II), and platinum(II) thiolate complexes, and darunavir (an antiretroviral drug used to treat and prevent HIV/AIDS) as a control drug (Figure 2B) to the inhibitor-binding domain of HIV protease was analyzed using UCSF Chimera (Pettersen et al., 2004;Kožíšek et al., 2014).While the 3D structure of HIV protease was retrieved from the Protein Data Bank [https://www.rcsb.org]with code 4LL3 (Kožíšek et al., 2014), the structures of the metal complexes were drawn using MarvinSketch (ChemAxon, 2013) and subsequently optimized and converted to 3D using Avogadro (Hanwell et al., 2012).The 3D structure of HIV protease was prepared by removing cocrystallized water molecules, non-polar hydrogens, and other non-standard residues using UCSF Chimera.The binding site for molecular docking was defined using a co-crystallized known inhibitor of HIV protease, darunavir (Figure 2).To map out the binding site, a grid box was used with dimensions centered at x = −8.58,y = 16.05, and z = 27.27 and size: x = 11.06,y = 6.35, and z = 16.01.To prepare for molecular docking, hydrogen and Gasteiger charges were added to the metal complexes, after which each compound was subsequently set as flexible, while the receptor was positioned in a rigid conformation using UCSF Chimera.The default exhaustiveness of 8 for AutoDock Vina was used to allow for a thorough conformation search and optimization.Compounds that showed the top three docking scores (mostly negative) were selected for further analysis via molecular dynamics (MD).

Molecular dynamics simulation of protein-ligand complexes
Molecular dynamics was performed using the GPU version of AMBER 18 with an incorporated PMEMD module (Le Grand et al., 2013) following standard protocols as elaborated in several computational studies (Ramharack et al., 2023;Agoni et al., 2018a;Agoni et al., 2018b;Olotu et al., 2019;Agoni et al., 2022).In all, four complexes generated from molecular docking were prepared for molecular dynamics simulation.The corresponding unbound (apo) form of the HIV protease was also prepared for simulation.
The compounds were parameterized using the ANTECHAMBER module, where atomic partial charges (AM1BCC) of Gaff, using the bcc charge scheme, were added (Salomon-Ferrer et al., 2013;Maier et al., 2015).For the gold atoms, parameters were derived using a combination of quantum mechanical calculations and established parameters from literature specific to gold (Mohammadnejad et al., 2015).Partial charges for the gold-containing ligands were assigned using the restrained electrostatic potential (RESP) method, following geometry optimization and electrostatic potential calculations at the HF/6-31G* level for non-metal atoms and a suitable level of theory for gold atoms (Schauperl et al. (2020).Bonding parameters for gold-ligand interactions were derived from quantum mechanical calculations and validated against experimental data when available.Nonbonding parameters, including van der Waals (vdW) interactions, were adjusted to accurately reflect the behavior of gold in a biological context.HIV protease was also parameterized using the Amber FF14SB force field (Salomon-Ferrer et al., 2013;Maier et al., 2015).Using the LEAP module, hydrogen atoms were also added, while the entire system was neutralized by adding counter ions (Na + and Cl − ), followed by a subsequent generation of ligand, protein, and complex topologies, as well as parameter files for each of the molecules.The systems were precisely solvated with water molecules using the TIP3P orthorhombic box with a size of 12 Å (Case et al., 2005).The solvated complexes were minimized initially for 2,000 minimization steps, applying a restraint potential of 500 kcal/mol, and then fully minimized for another 1,000 steps of the steepest descent without restraint.This was followed by the gradual heating of the systems from 0 K to 300 K for 50 ps, after which they were equilibrated for 500 ps, while the temperature and pressure were kept constant at 300 K and 1 bar, respectively (Berendsen et al., 1984).This was followed by MD production runs of 300 ns for each system, during which the SHAKE algorithm (Kräutler et al., 2001) was used to constrain all atomic hydrogen bonds.The MD simulation was initiated using a time step of 1 fs and coordinates saved at 1-ps intervals, followed by the subsequent analysis of trajectories using the integrated PTRJ and CPPTRAJ module (Roe and Cheatham III, 2013).The complexes and data plots were visualized using the graphical user interface of UCSF Chimera and Microcal Origin analytical software (Seifert, 2014).

Protonation state assignment
The accurate determination of protonation states for key residues in the active site is crucial for the fidelity of MD simulations and the design of effective protease inhibitors.HIV-1 protease contains two critical aspartate residues (Asp25 and Asp125) in its active site, which play a vital role in its catalytic mechanism (Adachi et al., 2009;Ghosh et al., 2016).Based on well-established practices and previous studies, these aspartate residues were modeled in their deprotonated forms in our simulations (Privat et al., 2020).This decision is supported by the typical pKa values of aspartate residues in enzymatic active sites, which suggest that they remain deprotonated under physiological pH conditions to act as nucleophiles in the catalytic process (Hofer et al., 2020).

Post-molecular dynamic simulation analysis
After a 300-ns MD simulation run, the coordinates of all the simulated systems were saved at 1-ps intervals, followed by the subsequent analysis of trajectories using the integrated PTRJ and CPPTRAJ module.The root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (RoG), and thermodynamic calculations were calculated and plotted for each of the resultant trajectories.The complexes and data plots were visualized using the graphical user interface of Discovery Studio (BIOVIA Discovery Studio, 2016) and Microcal Origin analytical software (Dassault Systèmes, 2016).

Thermodynamic calculations
The molecular mechanics/generalized born surface area (MM/ GBSA) method was used to calculate the binding free energy of the protein-ligand complexes (Miller et al., 2012;Genheden and Ryde, 2015).The binding free energy (ΔG bind ) of the protein-ligand complex was calculated as follows (Eqs 1-5): where ΔG bind is the summation of the gas phase and solvation energy terms minus the entropy (TΔS) term.
ΔE gas is the sum of the Amber force field internal energy terms ΔE int (bond, angle, and torsion), the covalent van der Waals (ΔE vdw ), and the non-bonded electrostatic energy component (ΔE elec ).The solvation energy is calculated as follows: G non polar γSASA + b. (5) The polar solvation contribution is represented as G GB , and G non polar is the non-polar solvation contribution and is calculated from the solvent-assessable surface area (SASA), obtained using a water probe radius of 1.4 A °.The surface tension constant (c) was set to 0.0072 kcal/mol, and b was set to 0 kcal/mol.
Although MM/GBSA methods are widely used for estimating binding free energies in molecular dynamics simulations, it is important to acknowledge their limitations, particularly when applied to metal-containing systems.One major limitation is the accuracy of solvation energy calculations for metal ions and metal-ligand complexes, which may not be adequately captured by continuum solvent models (Rastelli et al., 2010;Kongsted et al., 2009).Several studies have highlighted the challenges in accurately describing the solvation behavior of metal ions and metal-containing compounds using MM/GBSA methods.For instance, the non-polar solvation term in MM/ GBSA calculations relies on the SASA, which may not appropriately account for the complex solvation environments of metal ions (Wang et al., 2020).Additionally, the treatment of metal-ligand interactions, including electrostatic and van der Waals interactions, may require specialized parameterization and consideration of metal-specific effects, which are not fully addressed in standard MM/GBSA approaches (Genheden and Ryde, 2015).
All the new products are soluble in more polar and protic organic solvents, show good stability at room temperature, and are stable in air.Their stability and good solubility make them even more appealing candidates for biological screening since the stability of such compounds in solution is a vital consideration for biological evaluation.The complexes were characterized using spectroscopic techniques and microanalysis.The 31 P{ 1 H} NMR spectra of complexes 3 and 4 contain magnetically inequivalent phosphorus atoms, resulting in two peaks in each case, which resonated at 37.57 and 23.46 ppm for complex 3 and 40.16 and 30.32 ppm for 4, compared with the precursors 1 and 2 at 25.95 and 34.01 ppm, respectively.This is attributed to the cis and trans analogs of the products.Acetylated complexes of 3 and 4 were also obtained from a reaction of 1 or 2 and the acetylated gluconamido ethyl thiol ligand and isolated as an off-white solid of 5 or 6.The 31 P{ 1 H} NMR spectra of complexes 5 and 6 displayed a single peak in each case, which resonated upfield at 23.81 and 33.83 ppm, respectively, compared to the precursors.The phenyl carbon atoms display well-resolved signals in their 13 C NMR spectra.
In order to study the structural activity of these compounds, modified congeners of 5 and 6 were synthesized by reacting 2gluconamido ethyl thiol with compounds 1 or 2 in the presence of Et 3 N to yield complexes 7 or 8, according to Scheme 2. The 1 H and 13 C{ 1 H} NMR spectra showed the expected resonances.For example, the 31 P{ 1 H} NMR spectra of complexes 7 and 8 were measured at 29.25 and 34.67 ppm, respectively, which resonated downfield after the conjugation of acetylated 2-gluconamido ethane thiol with compounds 1 and 2. The micro-analysis results of complexes 3-8 were in agreement with the calculated values, which confirmed the complexation and purity of these complexes.

Molecular docking
To further establish the therapeutic potential of the mono and bimetallic glyco diphenylphosphino gold(I), palladium (II), and platinum (II) thiolate compounds synthesized, we explored their binding potential against HIV protease using molecular modeling techniques.This was premised on the previously reported inhibitory potential of gold(I) phosphine compounds against HIV-1 RT and HIV-1 PR by Fonteh and Meyer (2009).Molecular docking was performed, allowing for establishing the best binding poses of the compounds, which were subsequently scored for the ranking.The ligand-binding site of HIV protease was defined using its cocrystallized, known inhibitor.Docking scores that corresponded with compounds with the best pose were recorded and are given in Table 1.As shown, compounds 1, 2, 3, 4, 9, 10, and 11 showed higher binding energy scores than the standard inhibitor darunavir, suggesting the potential favorable inhibition of the enzyme.The top three potential inhibitors (compounds 2, 3, and 4) were selected for further modeling to provide detailed insights into the potential binding of these compounds to HIV protease.

Elucidating the structural binding mechanism of metal complexes relative to darunavir
The binding of darunavir to HIV protease for its inhibition is characterized by strong interactions between the 3(R),3a(S),6a(R)bis-THF moiety of darunavir and the main-chain atoms of aspartates 29 and 30 (King et al., 2004;Ghosh et al., 2006).These interactions account for the more favorable binding enthalpy of darunavir to HIV protease than that of amprenavir, a chemically related inhibitor broadly used in antiviral therapy.Darunavir was therefore used as a control to assess the potential binding mechanism of compounds 2, 3, and 4 toward HIV protease.Molecular docking allowed for the prediction of the binding modes of compounds within the inhibitor-binding site of HIV protease while also revealing the specific interactions that could inform their possible binding mechanism (Prieto-Martínez et al., 2018;Lin et al., 2020).

SCHEME 3
Synthesis of bimetallic complexes 9-11.As shown in Table 1, docking scores of −9.1, −9.4, and −9.9 kcal/ mol were calculated for compounds 2, 3, and 4, respectively.As a control, darunavir was redocked into its binding pocket, showing a docking score of −7.8 kcal/mol.Redocking of darunavir also allowed for an assessment of the reliability and validation of the docking protocol used.The docking scores calculated represented the bestbound conformations and provided insights into the likely binding affinities of the compounds (Dar and Mir, 2017;Prieto-Martínez et al., 2018).The compounds exhibited relatively higher docking scores than darunavir.The strong interactions of darunavir with ASP29 and ASP30, as previously reported, were highlighted, characterized by conventional hydrogen bonds and van der Waals interactions (Figure 3).The binding of darunavir also included conventional hydrogen interactions with ASP25 and GLY27.The binding of compound 4 was characterized by conventional hydrogen bonds with ASP25, ALA28, and ILE50 in several π-alkyl and alkyl interactions with other binding pocket residues, notably pi-alkyl interactions between CL and LEU23.Although the gold metal (Au) interacted with ARG8, it was unfavorable.Collectively, these interactions stabilize compound 4 within the binding pocket, resulting in a strong binding affinity.Compound 2 was also shown to interact with ASP25, ALA28, and ARG8, forming conventional hydrogen bonds in addition to π-alkyl and alkyl interactions with other binding pocket residues.Similarly, Au was shown to engage in an unfavorable interaction with ARG8.The binding of compound 2, on the other hand, involved van der Waals, π-alkyl, and alkyl interactions with binding pocket residues.The crucial residues ASP29 and ASP30 were shown to engage in van der Waals interactions with compound 2, whereas Au did not engage in any visible interaction, which accounts for lower binding affinities than those of other compounds.In general, the three primary metal complexes (compounds 2, 3, and 4) displayed interactions with comparable residues to darunavir, along with additional robust interactions with adjacent binding pocket residues.These supplementary interactions have the potential to bolster their binding affinity, as indicated by the docking calculation.

Exploring the conformational dynamics of HIV protease upon the binding of metal complexes
To elucidate the possible binding mechanism of the compounds toward the HIV protease, we assessed the impact of the binding of the top docked compounds (2, 3, and 4) by calculating the RMSD and RMSF of C-α atoms of HIV protease atoms over a 300-ns MD simulation.Estimating the RMSD and RMSF of individual atoms unravels the conformational or structural changes that could occur upon ligand binding.The therapeutic binding of chemical compounds to biological targets is typically associated with conformational and structural changes in both the ligand and target receptor, which influence the functions of the biological targets (King et al., 2004).As shown in Figure 4, the binding of all the compounds induced a reduction in the average RMSD and RMSF of the bound HIV protease relative to the unbound structures, suggesting that the binding of the compounds imposed some level of stability and decreased the residue flexibility of HIV protease.HIV protease was also shown to assume a more stable and rigid conformation when bound to the metal complexes than when bound to darunavir.The observed stability and decreased flexibility suggest that the binding of the metal complexes imposes structural rigidity, which could consequently impede the known flexibility of HIV protease, particularly its flap regions, a mechanism consistent with the binding mechanism of some known HIV protease inhibitors, such as darunavir.This observed structural dynamics could also consequently influence the binding interaction with active site residues and overall binding affinity as a relatively stable and rigid conformation could ensure the stability of crucial interactions required for the therapeutic modulation of HIV protease.Nonetheless, further studies are recommended to further validate these findings.Bound conformations and corresponding protein-ligand interaction dynamics of metal complexes, compounds 2, 3, and 4, and the known HIV protease inhibitor, darunavir.

Synthesized metal complexes exhibit favorable binding affinity toward HIV protease
We estimated the binding affinity of the top three synthesized metal complexes to HIV protease using the MM/GBSA-based binding free energy calculation approach (Ghosh et al., 2006).The MM/GBSA method considers several energy contributions, including vdW and electrostatic interactions, polar solvation (PS) energy, and non-polar solvent-accessible surface area energy.The degree of the binding affinity reflects the strength of interactions between the metal complexes and HIV protease and, making it a crucial component of their binding mechanism.The MM/GBSA calculations results, as given in Table 2, returned an estimated −18.81, −20.31, −26.73, and −31.17 kcal/mol for darunavir and compounds 2, 3, and 4, respectively.This result is consistent with the findings from our RMSD and RMSF calculations, in which the stable and rigid conformation observed in the metal complex-bound structures could have influenced the formation of relatively stable interactions to enhance binding affinities.The MM/ GBSA calculations provide additional molecular insights into the therapeutic potential of the synthesized metal complexes as they show enhanced structural modification of HIV protease.

Conclusion
A series of new gold(I) thiolate derivatives and three bimetallic complexes with amino phosphines and thiocarbohydrates as auxiliary ligands were synthesized and characterized.The structures and purity of complexes 1-11 were determined by NMR spectrometry and micro-analysis.It is believed that these thiocarbohydrate gold(I) complexes may show potential as anti-HIV agents since the composition of metals [gold(I), palladium(II), and platinum(II)], phosphines, and thiocarbohydrates is critical in maintaining the activity of the anti-HIV motifs.Molecular modeling further identified three of the metal complexes as potential HIV protease inhibitors, whose binding was characterized by strong binding affinity interactions with binding pocket residues, and the inhibition of the flexibility of the flap regions of the HIV protease was similar to that of the known HIV protease inhibitor, darunavir.This observation is based on (A) Comparative RMSD plots of unbound HIV protease (apo) and its complexed conformation when bound to compounds 2, 3, and 4 and darunavir.BSA with or without PVP over a 100-ns simulation.(B) Comparative RMSF plots of unbound HIV protease (apo) and its complexed conformation when bound to compounds 2, 3, and 4 and darunavir.BSA with or without PVP over a 100-ns simulation.The table shows average RMSD and RMSF estimations of simulated models over 300 ns.Frontiers in Chemistry frontiersin.org11 computational simulations using molecular dynamics, which suggest a potential mechanism of action for these compounds.However, it is important to note that these findings are theoretical and have not yet been confirmed through experimental validation.Although our computational results are promising, we did not perform experimental assays to directly demonstrate that the synthesized compounds act as protease inhibitors.Future work will focus on conducting biochemical assays, such as enzyme inhibition assays, to empirically validate the inhibitory activity of these compounds against HIV protease.

TABLE 1
Docking scores of the synthesized metal complexes toward HIV protease.

TABLE 2 MMPBSA
-based binding free energies of the synthesized metal complexes.ΔE ele , electrostatic energy; ΔE vdW , van der Waals energy; ΔG bind , total binding free energy; ΔG sol , solvation free energy; ΔG, gas phase free energy.